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Abstract. The electrostatic potential of any test charge distribution in Schwarzschild 
space with boundary values is derived. We calculate the Green's function, generalize 
the second Green's identity for p-forms and find the general solution. Boundary value 
problems are solved. With a multipole expansion the asymptotic property for the 
field of any charge distribution is derived. It is shown that one produces a Reissner- 
Nordstrom black hole if one lowers a test charge distribution slowly toward the horizon. 
The symmetry of the distribution is not important. All the multipole moments fade 
away except the monopole. A calculation of the gravitationally induced electrostatic 
self-force on a pointlike test charge distribution held stationary outside the black hole 
is presented. 
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1. Introduction 

In 1968, Israel proved the following theorem jl]: the Reissner-Nordstrom solution is the 
only static, asymptotically flat, electrovac solution of Einstein's equations for which the 
surfaces goo = constant are closed and simply connected and the event horizon g 00 = 
is regular. However, what happens, if a test charge is located near a Schwarzschild 
black hole? The electrostatic field of a point charge at rest in Schwarzschild space was 
derived by Hanni and Ruffini |2| and Cohen and Wald |§ for multipole fields and by 
Linet || in algebraic form. As the charge was slowly lowered into the hole, the authors 
found that the electric field of the charge remains well behaved, while all the multipole 
moments, except the monopole, fade away. From Israel's theorem they concluded that 
a Reissner-Nordstrom black hole was produced. Then in the early 1970s, when the 
black- hole uniqueness theorems were proved (particularly the theorems of Carter H, 
Hawking |J and Robinson 0), it became clear that an isolated black hole cannot have 
an electromagnetic field unless it is endowed with a net electric charge. 

The closely related problem of a magnetic field outside a compact magnetic star 
with a surface current was considered for dipole fields by Ginzburg and Ozernoi ||, and 
for multipole fields by Anderson and Cohen ||. Petterson [|Kj presented a calculation 



of the quasistatic axisymmetric magnetic field in a Schwarzschild background at radii 
both inside and outside the radius of the source. The authors found that the magnetic 
field vanishes for an observer at infinity when the source approaches the horizon. This 
is in accordance with a theorem by Price |TT|, which states that during the process of 



gravitational collapse, all electromagnetic multipole moments of the collapsing matter 
must disappear, except the electric monopole moment. However, the electric and 
magnetic fields become very large when the sources are very near the horizon. This 
fact is probably only a manifestation of the breakdown of the quasistatic treatment, 



and not a real effect. A time- dependent treatment by de la Cruz et all |12[ supports 
this view. They let a spherical shell of matter collapse at the speed of light upon a 
fixed magnetic dipole. The numerical result is that the magnetic field outside the shell 



decays to zero during the collapse. Wald [ 13 1 found the same sort of behavior when 



he calculated the electromagnetic field of an electrostatic or magnetostatic multipole 
of fixed strength placed at the centre of a massive, non-rotating, spherical shell. As 
the shell approaches its own Schwarzschild radius, all electrostatic and magnetostatic 
multipoles except the monopole decay to zero and the electromagnetic field remains 
finite on the shell. 

Electric and magnetic fields in the vicinity of massive compact objects are of 
particular interest for the study of black hole magnetospheres which give hope to 
understand the mechanism behind the enormous power outputs of active galactic nuclei 
and double radio sources. I shall not attempt to review the numerous works that 



have been written addressing black hole magnetospheres (see, e.g., |1J] and references 
therein). 

In section |2| we calculate the electrostatic potential of any test charge distribution 
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in a Schwarzschild background. For that, first, we calculate the Green's function for 
a spherical shell bounded by r = a and b and then, we make use of the generalized 
second Green's identity for p-forms (see the appendix). With this result, we are able 
to solve boundary value problems. An example is demonstrated in section |3|. As in the 
usual electrodynamics, we make in section [| a multipole expansion. It shows that for an 
observer far away from the sources all the multipole moments, except the monopole, go 
to zero as the charge approaches the horizon. This was known only for a point charge. 
The asymptotic properties, which are valid for any charge distribution, are presented 
in section g. A force is necessary in order to hold at rest a test charge distribution. In 
section ^ I recapitulate the works on the force and I show how one can recover their 
results with the Green's funktion method. There is a repulsive self-force and this force 
has for all static pointlike test distributions the same form. It depends only on the 
total charge Q. The derivation of the generalized second Green's identity for p-forms 
on (pseudo-) Riemannian manifolds can be found in the appendix. 



2. Electrostatic potential of a test charge distribution 

We use the standard Schwarzschild coordinates with the metric 

ds 2 = (l - ^ di 2 - (l - ^ 1 dr 2 - r 2 (dtf 2 + sin 2 tf dp 2 ) 



We write Maxwell's equations with differential forms [|15 



dF = , 5F = 4?r J . (2) 

The first equation in (^|) implies, by Poincare's lemma, the (local) existence of a potential 
A = A^dx^ 1 with F = dA. Then, the inhomogeneous Maxwell equations become 

5dA = An J . (3) 

Since the Maxwell equations are invariant under a gauge transformation of the potential 
A, we may make use of this gauge freedom and require the Lorentz condition 5 A = 0. 
Thus, Maxwell's equation (^|) may be written in the form 

DA = 4vr J , (4) 

where D:=(5od + do(5is the Laplace-Beltrami operator. In the static case, when 
the functions A^ are independent on t, the electrostatic potential A t decouples from 
Ai (i = r, <f). Since we are interested in the field of a static test charge distribution, 
we set the spacelike components of the current equal zero ji — (i — r, tp) . Then, we 
may take Ai = and obtain as the only non-trivial equation (/i = t) 

( 2m\ 11 1 

1 -J d r [r 2 A t , r ] + , % [ sin ■& A t ,t } + 2 . 2 - A t)VV = -4nj t ■ (5) 

\ r J r l L J rsini; r l sin v 

The comma in equation (M) denotes an ordinary derivative. 



Electrostatic boundary value problems in the Schwarzschild background 



As in the usual electrostatics, we need a Green's function in order to derive the 
general solution for any test charge distribution. Therefore, we consider two new 1- 
forms G = Gt dt and Sd = St dt with 



□G = 4tt(5 



D 



(6) 



We call G the 'Green form' and Sp the 'Dirac form'. St shall be the 5-function with the 
normalization J *Sd = 1- Thus 



■y 2m 



S{r - r') 6(& - #') % - ip') . 



(7) 



r 2 sin 

First, we solve the homogeneous equation of (|6]). If we make the separation ansatz 

G t (x, x') = Mr, r')Y lm ($, <p) , x= (r, 0, <p) , (8) 



l.m 



we obtain the following equation for Ri(r, r'): 



1 - 



2m 



d 
dr 



d^ 
dr 



-l(l + l)Ri{r,r') = 0. 



(9) 



The solutions of equation @ have been obtained independently bei Israel ]TJ and by 
Anderson and Cohen [[J. With the transformation u — — — 1 and the new function 
Ui(u) = a/(1 + u)/(l — u)Ri(u) equation 

1 



(1 -u 2 ) U[' (u) - 2uU[{u) + 



i becomes 
1(1 + 1)- 



U, u) = . 



(10) 



The solutions of dlOf) are the associated Legendre functions of the first and second kind 
Pi(u) and Q\(u) [jlB] . In the following, we take the way of Cohen and Wald 0. They 
give as the two linearly independent solutions of equation (|) 

for I = (Definition) , 



9i{r) 



Mr) 



2 l l\(l- l)\m l , ,dPi ( r 
W (r - 2m) dr"U 



for / ^ . 



(2/ + 1) 



-(r — 2m) 



dQi 



1 



'11a) 



[lib) 



2 l (I + l)!Z!m' +1 v ""' dr \m 

where Pi and Q\ are the two types of Legendre functions Then, they note three 
properties of gi(r) and fi(r), that will be important for the following analysis: 

(I) For I = 0, go(r) = 1 (by definition) and /o( r ) = l/ r - 

(II) For all /, as r — > oo, the leading term of gi(r) is r l , while the leading term of 
fi(r) = l/r l+l . 

(Ill) As r -> 2m, fi(r) — > finite constant, but d/i/dr blows up as ln(l — 2mr~ l ) for 
/ 7^ 0. Since gz(r) = (r — 2m) x (polynomial in r), so, as r — > 2m, gi(r) — > as 
(r — 2m) for / 7^ 0. 
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' i A l™ 9i(r) + B lm fi(r)] Y lm (#, <p) for r > r' 



Lm 



G t (x,x')={- (12) 
i C im 9i(r) + Am fi(r)} Y lm (#, <p) for r < r' . 

l,m 

The constants Ai m , Bi m , Ci m , Am can be determined by the boundary conditions. 
Suppose that the boundary surfaces are concentric spheres at r = a and b. The vanishing 
of Gt{x, x') for x on the surface implies 

Ci m gi{a) + Am fi{a) = for r = a , (13a) 
A lm g l (b) + B lm f l (b) = for r = b. (13b) 

-Ci m /fi(a) = D lm /gi(a) and b im := -A lm / f t (b) = B im / gi(b) equation (|TJ) 



With d lfn 
becomes 

G t (x,x') = I 



' bl ™ W)fi{r) - fi{b)gi{r)] Y lm (#, for r > r' 



l,m 



Y d im [gi(a)fi(r) - fi(a)gi(r)] Y lm (&, ip) for 



(14) 



r < r . 



Lm 



By continuity of at r = r' we obtain 

6bn [<?*(W) - fi(b)9i(r')} = d lm [a(a)/,(rO - /j(o)y,(rO] 

Let 

^Zm ^Zm 

: Pi(a)/,(r') - /,(a)^(r') = gW)W) ~ WoW) ' 



(15) 



ai;. 



then we have 



G t (£c, x 1 ) = 22 a imRi(r, r')Yi m ($, ip) 



Lm 



where 



Ri(r,r 



, N f [<7z(«) /z(r') - /z(a)<7z(r')][<?zWz(r) - /z(%(r)] r > ^ 



1 [9i(b)fi(r') - fi(b) gi (r')][ gi (a)fi(r) - fi(a) 9l (r)] for r < r> . 
We set equations fliTf), ( |T7| ) and (0) into (^) and multiply both sides by r 2 . Thus, 



(16) 



(17) 



Z,m 



2m \ d / 2 di2, 



r J dr \ dr 

47T 



-Z(Z + l)i? z (r,r') 



Y lm {$,y) 



sin?? 



(If 



In order to evaluate a; m , we multiply fli~8l) by sin?? YJ? m ,($, <p) and integrate over $ and 
Using the orthogonality condition of the spherical harmonics Y/ m (0, ip), we obtain 



«7r 



d / 2 dRA 1(1 + 1) 



dr \ dr 



i ^m 



2m" 



-47r5(r-r') W,^')- 



(19) 
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Integration over an infinitesimal interval at r' gives 



AirY l ^, < p') = a l , 



/2 



dRi(r > r') 



dr 



dRAr < r') 



— r 



,2 uni 



dr 



9i(r)-^:(r)- fi(r) — (r) 



= ai m r' 2 [gi(a)fi(b) - gi(b)fi(a)] 

= ai m r ,2 \jgi(a)fi(b) - gi{b)fi{a)} W{g h f h r') , 
where W{g\, fi, r') is the Wronskian of gi and fi at r'. The Wronskian 

Ui ... u n 



(20) 



W{r) 



u 



u 



(n-1) 



II. 



(n-1) 



(21) 



of n linearly independent solutions Ux(r),... ,u n (r) of a linear differential equation of 
the form u^' + a ri _ 1 (r)M < ^ n ^ 1 ^ + . . . + ao(r)w = 6(r) satisfies [[H] 

(22) 



(23) 



W(r) = W(r ) exp ( - / a n _i(s) ds 
In our case a„_i = ai(r) = 2/r and therefore 

2 



So, we have 

r 2 W(gi, fi, r) = constant . (24) 

We find the constant by evaluating W(gi, fi, r) for large values of r. Property (II) gives 
the values for the functions g\ and // for large r. Then, by equation (0), we have for 



all r 



-r\l + l) 



r 2 W(g h f h r) = r 2 

= - (2/ + 1) . 
Setting (p5|) in (|20|) and solving for ai m gives 



z-i 



Air 



1 



2Z + 1 </i(o)/,(6)-(/,(6)/i(a) 
Thus, the solution of @ is 



(25) 
(26) 



Lm 



(21 + 1) [1 - (9l(o)/l(6))/(/l(o)9l(6))] 



x ( 9l(r<) _|W /l(r<) )( /l(r>) _^| 9 , (r>) ) , (27) 

where r< (r>) is the smaller (larger) of r and r'. Now, we use the second Green's identity 
(equation (|A.17|) in the appendix) and set u := A(x), v := G(x,x') with A = A t dt and 
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G = G t dt. Since G t (x, x') vanishes on the surface dD, we obtain as the solution of (|5|) 

r t' 2 si n 
A t (x) = / jt(x') G t (x, x') g^r dr'd^'d^' 

J D 1 — 

_ J_ / dGt ^ x '^ r ' 2 sin tf'dtfV . (28) 

Since the functions are not dependent on time, we neglect the integration over t. One 
have to choose jt(x) in (p8|) in such a way, that 



t 2 sin ?9 

^drdtfd^ = Q, (29) 

D 1 — 

where Q is the total charge. 

In 1976, Linet presented an algebraic solution for a point charge situated at the 
point (t~o, $o, Vo) with r > 2m ||. He obtained his solution by modifying a particular 



analytic solution found by Copson Hl8| , which did not satisfy the boundary condition. 
Linet's solution is 

Alt \ ^ 

Af(x) = 

e (r - m)(r - m) - m 2 A(#, <p) 

r o r [(r — m) 2 + (r — m) 2 — m 2 — 2(r — m)(r — m)A($, </?) + m 2 A 2 ($, y?)] 1 ^ 2 
with 

A(i?, v?) = cos i? cos $o + sin •& sin $ cos(f — ipo) . (31) 

Now, it is easy to generalize equation (|30|) for any charge distribution j t (x). With (|A.17| ) 
again 

Aft*) = / *(*0 (a, *0 ~TTT dr'dW^' , (32) 



D 



1 r' 



1 (r — m)(r' — m) — 77i 2 A(i?, tp) 

r ' r [( r — m ) 2 + ( r ' — m ) 2 ~~ m2 ~ 2( r ~~ m){r' — m)A(i?, <£>) + m 2 A 2 (i9, y?)] 1 ^ 2 

A($, </?) = cos cos i?' + sin $ sin cos(y? — ip') . (34) 

The boundary values are equal zero. In order to consider some boundary values, one 
should add to (|33|) a solution of the homogeneous equation, which is fully determined 
by the choice of the boundary conditions. 



3. Solution of a boundary value problem 

As an example let us consider a concentric ring of charge of radius R > 2m and total 
charge Q located at $ — 1 7r inside a hollow sphere of radius b > R with the potential 
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Vi^d' , if 1 ) on the surface. With a — > 2m in (|27j), the derivative of G t , evaluated at r' = 6, 

flip) d 9l (r>) 



is 



9r' 



— 9i{r) 



l.m 



(21 + 1) 



dr' 



^,(6) dr' 



(35) 



Consequently, the potential of the hollow sphere inside r = b is, according to (|28 

V{#, <p')Y?J$\ tp') sin 



£ 

6 2 ^(r) / d/j(r' 



6 ^( & ) dr ' 



2/ + 1 V dr' 

The charge density of the ring can be read off from fl29] 
Q 1 - 2m 



r'=b. 



8{r - R) 8{$ - \k) . 



Y lm (&, <f) . (36) 



(37) 



2tt r 2 sin $ v 2 
Since the problem is axially symmetric, we set m = 0. Then, using fl27|) with a — > 2m 
equation (|28|) becomes 



n=0 

where we have used 



2 2 ™(n!) s 



(2n)\ 



2 2n (n\y 



, / \ /2n(^) / \ 

/an(r>) 77T92n{r > ) 

92n(0) 



for I = 2n , 
for Z = 2n + 1 . 



P 2n (cosi?) + $(a:) ,(38) 



(39) 



r< (r>) is the smaller (larger) of r and i?. 
4. Multipole expansion 

A located distribution of charge jt is non- vanishing only inside a sphere of radius R > 2m 
around the origin and the charge is in the space between r = 2m and R. Since we are 
interested in the moments outside the charge distribution, r < = r' and r> = r. With 
a — > 2m, b — > oo in (|27 ) and the vanishing of the potential on the boundary surface in 



(f2~8|) we find the following multipole expansion for the potential A t : 

A * = An E 2zTT /i(r) gL y5m( ^' ^ ' 

Z,77l 

where the multipole moments qi m are defined by 

/t' 2 sin i9' 
0,(0 *(r' l ^0 i _ g m J r/ dr' dtf' d^' . 

We give only the moments with m ^ 0, since the moments with m < are related 
through 

= (-l) m <?L • (42) 



(40) 



(41) 



Qi- 
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For I = and I = 1 equation ([II]) becomes 

<?oo = 4= > ( 43 ) 

V47T 

V Svr (44) 



where 




P= z-^xjtix)- — -^drd$d<p (45) 

J r 1 — 2m jr 

is the electric dipole moment and x = (x,y,z) with x = r sin cosy?, y = r sin^sin^, 

z = r cos$. For / = 2 we have 

1 / 15 

922 = 12 V 2tt ^ n + 2i ^ 12 ~~ ^ ' 

1 / 1 K 

to = -WiT (Q13 + 2Q23), (46) 

3 V o7T 

g2 ° = ^\/5 g33 ' 

where Qr,- is the traceless quadrupole moment tensor 

G« = / ( 3a ^ - V 2 ) *(*) i_2m/r ^ ^ ' (47) 
With equations (|4"3]) , (|44| ) and ( fl6D the expansion of i4 t (ac) P0| ) begins 

a^) = 5 + IM p . x + l J-^ Y, Q«*w + ■ ■ ( 48 ) 

Now, we go back to equation (|40"D . The electric field components for a given 
multipole can be expressed most easily in terms of spherical coordinates. In the local 
orthonormal basis 



/ 2m 1 

9° = Jl dt, 6 l = . = dr , 6 2 = rd$, 6 3 =rsm$dtp, (49) 

V r 2m/r 

the components of the field tensor F = dA with A = A t dt are 

F = -A t>r 6° A 9 1 - A t , d - 1 A 6 2 - A tip - 1 — L_ 0° A 9 3 , (50) 

- 2m/r r ^/l - 2m/r rsmtf 

-^01 = -^r , F 02 = -Etf , -Fo3 = • (51) 

The electric field with definite l,m has spherical components 

4vr /j(r) 9 

^ = " 2TTT gL rv /l ~ 2m/r ^ ' (52) 

771 4vr „ /,(r) im 

£ * = - 27TI *» 77i - 2, n/r SE? F "» w ■ 
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For a dipole p along the z axis, the fields in (|52| ) reduce to the form 

E r = -p ^ 1 }^ cos?9, E# = p — / l(y ^ sin#, ^ = 0. (53) 
dr 1 — 2m/r 

5. Asymptotic fields for r' — > 2m 

Cohen and Wald showed that, if one lowers a point test charge slowly toward r' = 2m, 
one produces a Reissner-Nordstrom black hole. Is this true for any charge distribution? 
First, we consider the case of a charge distribution at the point r' > 2m and an observer 
at 2m ^ r < r'. Thus, r< = r and r> = r'. Setting equations ([27]) and (|28|) (with 
a — > 2m, b — > oo) in (|50| ) the components of the field in the orthonormal basis (^) 
become 



'01 



02 



1 d#(r) 

~ 47r 2^ ^y— t Fim( ^ v) ^ 



2^ + 1 dr 



^2TTT^ (r) — M — /iim ' ( 54 ) 



•yl — 2m/r f 
r sin 



4vr ^ 1 01^(0,^) 

,171 



where 



/v* 2 si n 7?' 
j t (r>, tf, V ') f^) Y{ m {V', <fi i _ dr' dtf' . (55) 

For r near 2m, we see, by property (III) of the functions gi, that Fqi remains finite and 
Fq 2 ~ F 03 ~ 0[(1 — 2m/r) 1 / 2 ]. A stationary observer at r' > r « 2m sees a radial 
electrostatic field. Since j t can always have a term (1 — 2m/r)/(r 2 sin'#) (cf fl29|)), the 
corresponding term in fl5"5| ) vanishes. So, we see, by property (III) of the functions fi, 
that the field components F Qi (i = 1, 2, 3) remain finite at r = 2m as r' — > 2m. 

In the case r > r', we make use of equations ( |52"D and (f4"T|). As seen above, the term 
(r' 2 sin ■#')/( 1 — 2m/ r') in ( |4*T| ) vanishes with the help of the corresponding term in j t . 
From the fact that, for all I ^ 0, <#(r') — > as r' — > 2m (property (III)), we conclude 



that all the multipole moments of qi m (|4"I"|) except the monopole go to zero as r' — ► 2m. 



Since goo — Q/\f^ and /b( r ) = r \ we have the result that for all r > 2m 

-EV = ^ , Etf = E v = , as r — > 2m . (56) 

Thus, although the charge distribution does not possess any symmetry, the electrostatic 
field approaches the spherically symmetric Reissner-Nordstrom value E r = Q/r 2 for 
r' — > 2m. 



6. Force on a charge distribution 



One may ask what force is necessary to hold at rest a test charge distribution at a point 
(ro, $0; fo) outside a Schwarzschild black hole in a freely falling local system. A number 
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of authors [fT9-2~2f have investigated the problem by assuming that the gravitational 
field was weak; thus they worked to leading order in the small quantity m/r. Smith 
and Will [23] presented an exact calculation for a charged test particle held stationary 
near a Schwarzschild black hole. Zel'nikov and Frolov considered the influence of the 
gravitational field of a charged black hole on the self-energy of a charged particle [p^j 
and found that there is a mass shift of the particle in the gravitational field and that 
the absolute value of the mass shift coincides with the absolute value of the shift for a 



uniformly accelerated electron |25 



Smith and Will p3| got the result that the force is given by two terms. The first 
term is just the negative of the gravitational field that the hole exerts on the test particle. 
The second term is the gravitationally induced self-force of the particle. The self-force 
is repulsive and has the magnitude (in Schwarzschild coordinates) 

f- = ^ ■ (57) 

Since the hole is uncharged and the self-force vanishes as m — > 0, we must assume 
that the effect is induced by the spacetime curvature. The gravitational field obviously 
modifies the elctrostatic self-interaction of the charged particle in such a way that the 
particle experiences a finite self-force. 

For our calculation with the test charge distribution we use what was called the 
'global method' by Smith and Will in their paper | 23| . If we displace the charge slowly 
by a distance 5xq toward the hole then, according to the freely falling system, an amount 
of work SW is done given in this system by 

bW = -F ext ■ 6x . (58) 

Because of the gravitational red-shift, the energy SE received by an observer at 
asymptotic infinity is 

6E = y/g 0O (r )6W. (59) 

However, conservation of energy forces this energy to coincide with the change in the 
asymptotically measured mass — 5M of the system. In the freely falling system the 
coordinates can be chosen to be locally flat; that means 

9^ = V^, 9^9^ = 0- (60) 

The transformation between the arbitrary coordinates and the locally flat ones is done 
by 

^ = ^A^A^. (61) 
Hence, we have 

5x« = A™p 5x p . (62) 

Then, the external force on the test distribution in the locally flat system is 

rr 1 SM 
F^ xt = —== -= . (63) 

V0oo( r o) Arv6a% 
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The change in mass SM between two nearby, non-rotating black-hole configurations 



was calculated by Carter |26[ 

5M=^5A-^-5j G° V=g d 3 x + J G^h^V^g^x , (64) 

where x and A are the surface gravity and the area of the black hole, respectively, 
G^ v is the Einstein tensor, and h^ v is the difference in the metric between the two 
configurations. The integrals are to be evaluated over the exterior of the black hole. In 
our case, we can set 5 A = and ignore the term involving h^ u . Thus, with Einstein's 
equations, we have 

5M = -dJ T° oV ^d 3 x . (65) 

Since the energy-momentum tensor T^ v of the system has a mechanical contribution 
and an electromagnetic contribution, we split the 'energy' integral above into two terms 

- J T° y/=g~ d 3 x = U mcch + U cm . (66) 

With (|65"D and (]66f) equation ( |63D becomes 

_ 1 S(U mech + Uem) 

ext — / 7 r A 7T r ,, l D 'J 

V#oo( r o) A*vOX% 

To evaluate the mechanical contribution, we must choose a particle model. The easiest 
one is the model of an ideal fluid 

r = ( ? °+j))«Y-|)f , (68) 
where p is the pressure, g° the density and the velocity field with 

^itV = 1 . (69) 
Since for the particle density g° 

g y/=g~ d 3 x = m , (70) 

where mo is the total mass of the charge density, we obtain 

Q° = ^L5 3 (x-x ). (71) 

In our system p — 0, thus 

T>»=J?**f*fp( x - Xo ), (72) 
yj—g dr dr 

where x M (r) is the worldline of the charge density and r the proper time. We write ( |T2]) 
in manifestly covariant form 

T ,u = m^fdx^dx^ ^ _ Xo(r))dr (73) 
yj-g J dr dr 

In the local flat system u M = (1, 0, 0, 0). One finds from that 

dr = v ^dt. (74) 
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Choosing r = when t = 0, we have 

1 



t, x ,y ,zo 



(75) 



Now, the mechanical contribution to the energy is easily found from (^) and (|75|) to be 

f^mech = -mO"\Atoo( r o) ■ (76) 

To evaluate the electromagnetic contribution, we first note that 



rpfj,U ^ 



1 



F n F \v + _ g y.u F <r\ Ftr 



4 



So, it follows that 



r 

L em o 



07T 



and hence 



U em = ±- f gMg^dtAodjAo^g-^x 

87T J 

Via an integration by parts, one obtains 



(77) 

(78) 
(79) 



~ [d 3 xA d l (^g-g 00 gVd J A ) + ±- [ d 2 S, 
8tt J 8vr 

OTT Jr=2m 



-gg 00 g^A,d 3 A 



(80) 



The surface integrals vanish because drops off as 1/r at large distances and the 
angular integration on the horizon averages to zero as one can see using (|3~0"D. One is 
left with the first term, which because of (fj) becomes 

1 r j°AoV=gd*x. (81) 



12) 



Now, we replace Aq with the help of our Green's funktion G^(x, x') (equation (R3 

1 



fix 



IV^ft / J°(x')^gJ^)d 3 x'Gf(x,x') 



In the integral above, Gf(x, x') corresponds to the energy integral of a point test charge 
at x = x'. We can this clarify by setting j°(x') = J d 3 x" —g{x") 5°(x' — x")j°(x"). 
Then, equation (82) becomes 



with 



ui 



J f(x)y^gjxjd 3 x J f(x")^gJ^)d 3 x"U s cm 



l - j d 3 x'^g\^)5\x' -x") Gt(x,x') . 



53) 



34) 



Gt(x, x') has a singularity at the point x = x' . In order to separate out the divergent 
behavior, we expand Gf(x, x') about x = x' . We set \x — x'\ = a and take the limit 
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a — > after the integration. This corresponds to giving the point a finite radius a. For 
the integration it proves convenient to write G^(x,x') in isotropic coordinates 

g = | [r — m + (r 2 — 2mr) 1 / 2 ] , 
x = £> sin cosy? , 
y = g sin $ sin (p , 

Z = pCOS'i? . 



i5) 



Then 



e 1 + 



m 
2g~ 



1/2 



and the Schwarzschild metric has the form 
g = h 2 (\x\)dt 2 -f 2 (\x\)dx 2 

with 

l-m/2g 
w l + m/2e' M ; 

After a little calculation G^(x,x') (equation (^)) becomes 



1 + 



m 



2<? 



(86) 
(87) 
(88) 



£'(l + m/2£') 2 £(i + m/2^) 5 



m + g 



, ( q 2 - ^( x ' x + y'y + z ' z ) + q' 2 



1/2 



g 2 — 2(x'x + y'y + z'z) + g' 2 

\ 1/2" 
12 x ' 



+ ~q' 



g 2 — 2(x'x + y'y + z'z) + g' 
£ 2 - 0s(x'x + y'y + z'z) + g 



n'2 



19) 



where g = m 2 /4:g. For x' = (0,0,6) we get equation (39) of [23|. Now the expansion 
about x = x' gives 

1 l-m/2g' 



Gf(x,x' 



x - x'\ (1 + m/2g') 3 

— g' + m(l — m/4g') 
+ g' 3 (l + m/2g')(l - m/2g') 
m 1 



g' 2 (1 + m/2g'Y 



1 1 - m/2g' 

1 ; X (X — X 



x'(x - x') + 0(x - x') 2 

O(x-x') 2 



(90) 



g' 2 1 + m/2g' 

In equation (^4|) is \J—g(x') 5° = 5 3 (x' — x"). Like for Gf(x, x'), we build a ball around 
the point x'. Both, x" and x, are in a neighbourhood of x'. So we can set x = x" in 
the expression of 5° and substitute 

\/-g{x') 5° = -— lim 6(\x - x'\ - a) . (91) 



Using the expansion 



Ui 



yields 



47T0 2 a^O 

the integral 
1 l-m/2g' 

2a (1 + m/2g'f + 2^' 2 (1 + m/2g') 4 



m 



(92) 
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Now, we have the result for U em (|83|), which is valid for any charge distribution. 

If we choose pointlike distributions such as point charges or a ring of charge (cf, 
section |3p, j° has the form 

V^gf = Q6 3 (x-x ) , 

where Q is the total charge. Then, we obtain from the equations ( |TBj ) and 

1 



(93) 



^mech U c 



l—m/2g Q 2 l — m/2g mQ 2 
-m — - h — t^— - — — + 



1 + m/2g ' 2a (1 + m/2g ) 3 ' 2g 2 (1 + m/2g )* ' ^ 

Now, in order to renormalize the mass m to M, we must express the radius of the ball 
charge in local freely falling coordinates. With x M = x u we obtain 

a = a(l +m/2g Q ) 2 . (95) 

Then, we can write 



mech 



+ U cm = M 1 - m/2go + mQ2 



1 + m/2g 



2eg 



i + m/2^) 4 



with 



M 



-mi + lim — . 

5-0 2a 



(96) 



(97) 



It follows from ( |67| ) that 

7 Mmxl 



mQ 2 x l c 



1 



ext 



;i + m/2^o) 3 1 - m/2 eo 



1 + m/2^> c 



(9* 



For #0 = (0, 0, b) equation (j98|) is in agreement with equation (75) of [[23 
Schwarzschild coordinates (|98l) becomes 

' sin $ cos fo \ 
sin sin ip 

cos^o 



Mm 



ext 



2m 

To 



-1/2 



mQ 2 



\ 



In 



(99) 



/ 



Equation (|98|) (and (j99|) ) holds only for pointlike distributions and if we assume for 
U mec h an ideal fluid model. However, the repulsive self-force depends not on the chosen 
mechanical model, but only on the chosen charge distribution. Thus, the electrostatic 
self-force for pointlike test distributions, i.e. the second term on the right-hand side of 
equation (|98|) or (|99|) , is always (in the freely falling system) 

FL = , (100) 



where F e r xt 



is the radial component of the self-force in terms of the pointlike test 
distribution's Schwarzschild radial coordinate. The force is radially directed because 
of the spherical symmetry of the spacetime and depends only on the total charge Q. 
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We generalize the wellknown Green's identities for p-forms on a (pseudo-) Riemannian 
manifold. Loomis and Sternberg P7| and Flanders [|2S| give the derivation for functions 
on the n-dimensional Euclidean space. The result for Riemannian manifold can be found 
in section 21 of Holmann and Rummler's book [29[|. Thirring also presents a version of 
Green's second identity [0O| . However, his result is not well adapted to our case. Since 
the Schwarzschild spacetime is pseudo-Riemannian, we have to generalize the formula 
in p9 



8D 



du A * dv + (-1) 



D 



u A d * dv . 



(A.l) 
(A.2) 

(A.3) 



D 



Let (M, g) be an n-dimensional oriented (pseudo-) Riemannian manifold and let D 
be a region of M with smooth boundary such that D is compact. For the two forms 
a G /\ p (M) and (3 G f\ q (M) the anti-Leibniz rule gives 

d(a A f3) = da A (3 + (-l) p a A d/3 . 
Now, let u, v G f\ p (M) be two p-forms. Then, with ( |A.1| ), we obtain 

d(u A * dv) = du A * dv + (—l) p u Ad* dv . 
Using Stokes theorem we obtain 

u A * dv = f d(u A * dv) 

Note that dwA* dv = dvA* du. If we write down ( |A.3|) again with u and v interchanged, 
and then subtract it from ( |A.3j ), we have 

/ (m A * - t; A * du) = (-l) p / (u A d * dv - v A d * du) . (A.4) 

J 3D J D 

Now, we write the right-hand side of (|A.4fj in a different form. For every form 
u G A fc (M) (see 0) 

* *cu = (-l) k{n - k) sgn(g) u . 

This gives us 

(_!)*(*-*) sgn (0) * * u = u 

With uj = d * dv and k = n — p 

d*dv = (_i)(™-p)(«-p-«) S g n (^) * * d * dv 

= (-I)p(p-") sgn(^) * * d * dv . (A.7) 
The codifferential S : /\ q (M) -> A g -i( M ) is defined by p| 

5 := sgn(^) (-1)^+™ * d * . (A.8) 

(A.9) 



(A.5) 
(A.6) 



We solve (|A.8| ) for *d* 

*d* = {-l)- nq - n sgn(g) 6 . 
Setting ( |A.9| ) in ( |A.7|) we obtain for d* dv (and in analogous way for d* du) (q — p + 1) 

d*dv = (-l) p *5dv, d*du = (-l) p *8du . (A.10) 
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(it A * dv — v A * du) = / (u A * 5dv — v A * 5ckx) . 
<9r> </d 

Since for two p-forms a, (3 

a A* (3 = (3 A* a , 
equation (|A.11|) becomes 

[u A * dv — v A * du) = / (5dv A * u — 5du A * v) . 
) Jd 

Now, we consider the combination 

5u A *v G j\ n -i(M) . 

Using the anti-Leibniz rule and Stokes theorem again we have 
5uA*v = I d5u A * v + (— 1) 

i 

We interchange u and v and subtract the new equation from ([A. 13 ) 



dD 



p ~ l / 5u Ad* v . 
dD Jd Jd 



(5u A *v — 5v A * u) = / (d5u A * v — d5v A * u) 
dD Jd 



+ (-ir 1 / 

Jd 



u A d * v — 5v Ad* u) 



Now (by definition (JO])) 

5u Ad* v = sgn(g) (_l) n P+ n * d* u A d* v 

= sgn(g) (— l) n P +n *d*vAd*u = 5vAd*u. 
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(A.11) 



(A.12) 



(A.13) 



(A- 14) 



(A.15) 



If we set ( |A.15| ) into ( |A.14| ), the two last terms on the right-hand side of ( |A.14j ) cancel 
and we find 



(5v A * u — 5u A *v ) = / (d5v A * u — d5u A * v ) . (A.16) 
IdD Jd 

Now, we add (|A.12|) and (|A.16|) and obtain Green's second identity written with p-forms 

on a (pseudo-) Riemannian manifold 



(5v A*u — 5uA*v + uA*dv — v A * du) = / (Ov A * u — □« A*v) 
'dD Jd 

where □:=do5 + 5odis the Laplace-Beltrami operator. 



(A.17) 
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